Compact and Flexible Basis Functions for Quantum Monte Carlo Calculations 
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Molecular calculations in quantum Monte Carlo frequently employ a mixed basis consisting of con- 
tracted and primitive Gaussian functions. While standard basis sets of varying size and accuracy 
are available in the literature, we demonstrate that reoptimizing the primitive function exponents 
within quantum Monte Carlo yields more compact basis sets for a given accuracy. Particularly large 
gains are achieved for highly excited states. For calculations requiring non-diverging pseudopoten- 
tials, we introduce Gauss-Slater basis functions that behave as Gaussians at short distances and 
Slaters at long distances. These basis functions further improve the energy and fluctuations of the 
local energy for a given basis size. Gains achieved by exponent optimization and Gauss-Slater basis 
use are exemplified by calculations for the ground state of carbon, the lowest lying excited states of 
carbon with 5 S° , :i P°, 1 D° , :i F° symmetries, carbon dimer, and naphthalene. Basis size reduction 
enables quantum Monte Carlo treatment of larger molecules at high accuracy. 



I. INTRODUCTION 

In traditional quantum chemistry (QC) calculations, 
molecular orbitals are often expanded in a combina- 
tion of contracted Gaussian basis functions and primi- 
tive Gaussian basis functions. For each occupied orbital, 
a contracted function is constructed to reproduce the 
corresponding atomic orbital from an effectively single- 
electron theory such as Hartree-Fock (HF) Jn, or the 
natural orbital from a post-HF method [!, [J 

While a single primitive Gaussian has incorrect long- 
range asymptotic behavior, a contracted basis function 
can reproduce the correct asymptotics over a reasonable 
range. However, even contracted functions are unable to 
produce the correct electron-nucleus cusps @ since they 
have zero gradient at the origin. Despite these shortcom- 
ings, Gaussians are used in QC calculations because they 
p_ermit analytical evaluation of the two-electron integrals 

In contrast to traditional QC methods, quantum 
Monte Carlo (QMC) calculations Q enjoy greater wave- 
function flexibility by using Monte Carlo integration to 
evaluate matrix elements. In particular, bases need not 
be restricted to Gaussians. For calculations employing a 
potential that diverges at the nucleus, Slater basis func- 
tions can exactly reproduce the correct electron-nucleus 
cusp and long-range asymptotic behavior of the orbitals. 
In fact, for all-electron QMC calculations, highly accu- 
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rate results have been obtained by employing compact 
basis sets consisting of Slater functions with optimized 
exponents 

Conversely, the basis sets used for non-divergent pseu- 
dopotential calculations in QMC have deviated little 
from typical QC basis sets. For these pseudopoten- 
tials, orbitals have no electron-nucleus cusp. In this 
case, Gaussian basis functions are more appropriate than 
Slater functions at small electron-nuclear distances but 
still have incorrect long-range asymptotics. 

Contracted and primitive Gaussian functions are fre- 
quently splined on a radial grid for QMC. Splining con- 
tracted Gaussians presents a definite computational ad- 
vantage since evaluating polynomials is much cheaper 
than evaluating a linear combination of Gaussians. In 
contrast, splining primitive Gaussians provides minimal 
benefit at best. 

We propose two ideas for improving basis sets for pseu- 
dopotential calculations in QMC. First, primitive basis 
function exponents are optimized for each system. This 
provides greater accuracy with a compact basis for a wide 
range of chemical environments and excitation levels. To 
facilitate optimization, the primitive basis functions re- 
main analytic while the contracted functions are splined. 

Second, we propose a novel form of primitive ba- 
sis function appropriate for calculations involving non- 
diverging pseudopotentials. These primitives, which we 
call Gauss-Slater (GS) functions, have the short-range 
behavior of a Gaussian function and the long-range be- 
havior of a Slater function. 

The utility of our improvements is demonstrated by 
calculations for carbon, the lowest lying excited states 
of carbon with 5 S°, 3 P°, 1 D°, 3 F° symmetries, carbon 
dimer, and naphthalene. 
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This paper is organized as follows. In Section UH the 
form and properties of Gauss-Slater functions are intro- 
duced. In Section IIII| results of our calculations are dis- 
cussed. In Section [IV! concluding remarks are provided. 
In the Appendices, technical details are discussed. 



II. GAUSS-SLATER BASIS FUNCTIONS 

We define Gauss-Slater (GS) functions as 

where r, 9, <p are the standard spherical coordinates, n is 
the principal quantum number, I is the azimuthal quan- 
tum number, m is the magnetic quantum number, is 
the normalization factor, and ZJ n (9, 4>) is a real spherical 
harmonic. 

Notice that for r <C 1 the GS behaves like a Gaussian: 

fiiJr, 9, 4>) - N< r^e-W zp(9, 0), (2) 

and for r ^ 1 the GS behaves like a Slater: 

vii m (X, 9, 4>) = Ni r n ~ l e-< r Z?{9, 4>). (3) 

The GS drift velocity and local energy are well behaved 
at long distances, while for Gaussians these quantities 
diverge. 

Unlike Gaussians and Slaters, normalization of GSs has 
no closed form expression. Nevertheless, normalizing an 
arbitrary GS is trivial with the following scaling relation 
(see Appendix lAl between N% and iV 1 , 

N C = C+ l/2 ^ (4) 

Since GSs are not analytically integrable, the expo- 
nential part must be expanded in Gaussians for use in 
quantum chemistry programs that employ analytic inte- 
grals for evaluating the matrix elements. This expansion 
is 

where is the i th expansion coefficient, and is the 
i th Gaussian exponent. Additionally, the following scal- 
ing relations (see Appendix [A| hold for the expansion 
coefficients and Gaussian exponents: 

<4 = cV (6) 

4 = c\. (7) 

Once the Gaussian expansions are found for unit expo- 
nents, expansions of arbitrary GSs follow immediately 
from the scaling relations. 



III. RESULTS 

For all applications discussed in this paper, variational 
Monte Carlo (VMC) and diffusion Monte Carlo (DMC) 
[HI calculations are performed with the CHAMP QMC 
code [l(J and employ the pseudopotentials and accom- 
panying basis sets of Burkatzki, Filippi and Dolg (BFD) 
PH . We choose these pseudopotentials and basis func- 
tions since they were constructed for use in QMC and 
have proved to be quite accurate. 

The wavefunction is of the standard Slater- Jastrow 
form. All wavefunction parameters including Jastrow 
parameters, Configuration State Function (CSF) coeffi- 
cients (where applicable), orbital coefficients, and primi- 
tive exponents (where applicable) are optimized via the 
linear method [12Hl4 |. Optimization is performed on a 
linear combination of the energy and variance of the local 
energy with weights 0.95 and 0.05, respectively. Opti- 
mizing just the energy yields slightly lower energies and 
somewhat higher variances. 

For each system considered, calculations are performed 
with three different basis sets: (1) the BFD basis, (2) 
fixed contracted functions and analytical Gaussian prim- 
itives with optimized exponents, and (3) fixed contracted 
functions and analytical Gauss-Slater primitives with op- 
timized exponents. We refer to these basis sets as BFD, 
G, and GS, respectively. 

These three cases allow us to evaluate the improve- 
ments our two methods provide to the current basis sets 
used in QMC. First, if both the G and GS basis sets 
significantly outperform the BFD basis, then the utility 
of reoptimizing the basis exponents within QMC will be 
established. Second, the utility of the GS basis depends 
on its performance relative to the G basis. 



A. Ground State Carbon Atom 

For the carbon atom ground state, 3 P, we consider a 
complete active space (CAS) wavefunction with an active 
space generated by distributing the four valence electrons 
among the thirteen orbitals of the n = 2 and n = 3 shells. 
Denoted by CAS (4, 13), this wavefunction consists of 83 
CSFs comprised of 422 determinants. 

In general, a single Slater determinant will not be a 
CSF when a certain number of electrons have been ex- 
cited relative to the ground state HF Slater determinant. 
However, a CSF can be produced from an arbitrary Slater 
determinant by applying projection operators for angu- 
lar momentum L and spin S. Since states with the same 
L and S but different L z and S z are degenerate, we are 
free to choose convenient L z and S z states. We choose 
L z = to make the wavefunctions real to within a posi- 
tion independent phase, and we choose S z = S to yield 
the minimum number of determinants in the CSF. Since 
the carbon ground state has L = 1, S = 1, the projection 
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operators are of the form 



(8) 
(9) 



where the product over all possible angular momentum 
and spin values omits the desired L = 1 and S = 1 values. 

Carbon atom VMC results for each basis set are shown 
in Table HI Included for comparison, coupled cluster cal- 
culations with single and double excitations and pertur- 
bative triple excitations (CCSD(T)) values for the BFD 
basis [HI exhibit much larger dependence on basis size 
than QMC results. 

Both the G and GS basis sets outperform the BFD ba- 
sis set. The 2z G basis exhibits a modest gain of 0.3 mH 
in energy compared to the corresponding BFD basis. The 
2z GS basis exhibits larger gains of 1 mH in energy and 
28 mH in a, the root-mean-square (RMS) fluctuations of 
the local energy. The 3z GS basis yields identical results, 
within statistical error, to the hz BFD basis. 

Carbon atom DMC results for each basis set are shown 
in Table [TTJ These calculations were performed with a 
time step of r = 0.01 H _1 which leads to negligible time 
step error for these high quality wavefunctions. DMC 
depends less on basis size than VMC, as is immediately 
apparent from the data. Nevertheless, both the G and 
GS basis sets outperform the BFD basis set. The 3z GS 
basis yields identical results, within statistical error, to 
the hz BFD basis, and an energy 0.1 mH lower than the 
3z BFD basis. 

Both VMC and DMC results indicate that reoptimiz- 
ing primitive basis function exponents provides improve- 
ments which can be significant for the GS basis. In large 
systems, the ability to use a 3z basis in place of a 4z or hz 
basis determines whether a calculation can be performed. 



TABLE I: VMC energy and RMS fluctuations of the local 
energy, a, in Hartrees for CAS (4, 13) ground state of carbon 
using BFD, G, and GS basis functions. Statistical errors on 
the last digit are shown in parentheses. For each n, the nz 
basis consists of n S functions, n P functions, and n — 1 D 
functions. CCSD(T) values for the BFD basis are included 
for comparison [111 ]. 



lype 


oize 


jjjIltJI {J y in ) 


a (ft) 


BFD 


2z 


-5.43161(3) 


0.1395(6) 




3z 


-5.43306(2) 


0.099(3) 




42 


-5.43332(2) 


0.0904(2) 




52 


-5.43341(2) 


0.0905(4) 


G 


22 


-5.43196(3) 


0.138(2) 




32 


-5.43324(2) 


0.0989(5) 


GS 


22 


-5.43264(2) 


0.1114(4) 




32 


-5.43344(2) 


0.0898(2) 


CCSD(T) 


22 


-5.409230 


N/A 




32 


-5.427351 


N/A 




42 


-5.431486 


N/A 




52 


-5.432494 


N/A 



TABLE II: DMC energy in Hartrees for CAS(4,13) ground 
state of carbon using BFD, G, and GS basis functions. Sta- 
tistical errors on the last digit are shown in parentheses. For 
each n, the nz basis consists of n S functions, n P functions, 
and n — 1 D functions. Calculations were performed with a 
time step of r = 0.01 H~ which leads to a negligible time 
step error for these high quality wavefunctions. 



Type 


Size 


Energy (H) 


BFD 


22 


-5.43314(2) 




32 


-5.43395(2) 




42 


-5.43404(1) 




52 


-5.43407(1) 


G 


22 


-5.43342(2) 




32 


-5.43400(2) 


GS 


22 


-5.43356(2) 




32 


-5.43407(1) 



B. Excited States of Carbon 

We consider the lowest lying excited states of carbon 
with 5 S°, 3 P°, 1 D°, and 3 F° symmetries. These states 
have configurations 2s 1 2p 3 , 2s 2 2p 1 3s 1 , 2s 2 2p 1 3d 1 , and 
2s 2 2p 1 3d 1 , respectively. The 3 P°, 1 D°, and 3 F° states 
have much higher energy than the ground state and 5 S° 
excited state. 

The dominant CSF for each of these three states oc- 
cupies orbitals that are unoccupied in the HF ground 
state. For fair comparison, the BFD basis therefore must 
be augmented. The diffuse functions of the aug-cc-pVnZ 
basis sets [4j, |lg, ll7j are used for this purpose. The BFD 
nz basis then becomes an (n + l)z basis. 

Obtained by application of the projection operators 
discussed in Section IIII Al the dominant CSF for each 
of the four excited states has one, one, four, and three 
Slater determinants, respectively. 



VMC results for energies and a of each system are 
shown in Figures [1] and [5] In all cases, the reoptimized 
exponents provide significant gains in both energy and 
a. Results for the three higher lying states demonstrate 
that reoptimized exponents are essential for describing 
states containing orbitals unoccupied in the HF ground 
state. For these systems, 2z results using the G and GS 
basis sets are substantially better than hz results for the 
augmented BFD basis set. In the most extreme case of 
3 F°, the 2z GS basis results in 30 mH lower energy and 
110 mH lower a than the hz BFD basis. 

The importance of the reoptimized exponents is evi- 
dent for the excited states of carbon. However, benefits 
of the GS basis relative to the G basis are never more 
than several tenths of a mH. On the scale of the plots in 
Figures [1] and O many G and GS basis results coincide. 
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FIG. 1: VMC energies in Hartrees for lowest lying excited 
states of carbon with 5 S°, 3 P°, 1 D°, 3 F° symmetries. For 
each n, the nz basis consists of n S functions, n P functions, 
and n — 1 D functions (where applicable). For 3 P°, 1 D°, 
S F° calculations, the BFD basis is augmented with diffuse 
functions of the aug-cc-pVnZ basis sets [J, [lf^ [TtJ . For G and 
GS basis sets, only 2z and 3z calculations were performed. In 
many cases, results for G and GS bases are indistinguishable 
on this scale. 



These calculations were performed with a time step of 
t = 0.005 H _1 which leads to negligible time step error. 
The 2z G and GS basis sets significantly outperform the 
corresponding BFD basis set. The 2z GS basis yields a 
result that is essentially converged with respect to basis 
size. 

Both VMC and DMC results indicate that reoptimiz- 
ing primitive basis function exponents provides improve- 
ments which can be significant for the GS basis. 



TABLE III: VMC energy and a in Hartrees for single deter- 
minant carbon dimer ground state. Statistical errors on the 
last digit are shown in parentheses. For each n, the nz basis 
consists of n S functions, n P functions, n — f D functions, 
n — 2 F functions, and n — 3 G functions. 



Type 


Size 


Energy (H) 


a (H) 


BFD 


22 


-11.02644(4) 


0.4343(9) 




32 


-11.03003(4) 


0.4172(4) 




42 


-11.03094(4) 


0.4127(7) 




52 


-11.03095(4) 


0.4113(6) 


G 


22 


-11.02707(4) 


0.4288(7) 




32 


-11.03030(4) 


0.4183(6) 


GS 


22 


-11.02968(4) 


0.4191(6) 




32 


-11.03065(4) 


0.4109(6) 
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FIG. 2: RMS fluctuations of VMC local energies. See Figure 
[T]for notation and details. 



TABLE IV: DMC energy in Hartrees for ground state of car- 
bon dimer using BFD, G, and GS basis functions. Statistical 
errors on the last digit are shown in parentheses. For each n, 
the nz basis consists of n S functions, n P functions, n — 1 
D functions, n — 2 F functions, and n — 3 G functions. Cal- 
culations were performed with a time step of r = 0.005 H _1 
which leads to negligible time step error. 



Type 


Size 


Energy (H) 


BFD 


22 


-11.05561(3) 




32 


-11.05719(4) 




42 


-11.05728(4) 




52 


-11.05723(4) 


G 


22 
32 


-11.05632(3) 
-11.05717(3) 


GS 


22 


-11.05702(4) 




32 


-11.05719(3) 



C. Carbon Dimer 

Single determinant calculations were performed for the 
carbon dimer with initial wavefunctions generated from 
the QC code GAMESS (3- 

VMC results for each basis set are shown in Table Mil 
The G and GS basis sets outperform the BFD basis set. 
In particular, the 2z G basis attains a 0.6 mH lower en- 
ergy than the corresponding BFD basis, and the GS basis 
yields a 3.2 mH lower energy than the BFD basis. The 
3z GS basis yields an energy within 0.3 mH of and a a 
identical to the hz BFD basis results. 

DMC results for each basis set are shown in Table HVl 



D. Naphthalene 

Single determinant calculations were performed for 
naphthalene, CioHs, with initial wavefunctions generated 
from the QC code GAMESS [3. 

Calculations were performed only for the 2z basis. The 
intention of this section is not to produce an energy con- 
verged with respect to basis size, but to demonstrate that 
the utility of reoptimizing primitive basis functions, and 
GS primitives in particular, extends to large systems. 
VMC and DMC results for each basis set are shown in 
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Table [V] The DMC calculations were performed with a 
time step of r = 0.01 H -1 . 

At both the VMC and DMC level, wavefunction qual- 
ity increases by reoptimizing the primitive Gaussian ex- 
ponents. The GS basis provides further improvement. In 
particular, even at the DMC level, the 2z G basis attains 
a 4 mH lower energy than the corresponding BFD basis, 
and the GS basis yields a 15 mH lower energy than the 
BFD basis. This is significant since DMC is less sensitive 
to basis set choice than VMC. 

TABLE V: VMC energy and a, and DMC energy in Hartrees 
for single determinant ground state naphthalene, CioHs. Sta- 
tistical errors on the last digit are shown in parentheses. DMC 
calculations were performed with a time step of r = 0.01 H _1 . 
For carbon, the 2z basis includes 2 S function, 2 P function, 1 
D function. For hydrogen, the 2z basis includes 2 S functions 
and 1 P function. 



Type I Size I VMC Energy (H) | VMC a (H) | DMC Energy (H) 
BFD 2z -61.5193(5) 0.980(1) -61.6479(5) 

G 2z -61.5273(4) 0.938(1) -61.6518(5) 

GS 2z -61.5438(4) 0.927(2) -61.6634(5) 



IV. CONCLUSION 

Quantum Monte Carlo (QMC) methods have the ad- 
vantage over standard quantum chemistry (QC) methods 
of rapid convergence with increasing basis size. Basis- 
size dependence for pseudopotential calculations in QMC 
is further reduced by two basis set improvements intro- 
duced in this work. Calculations for the ground state 
of carbon, the lowest lying 5 S°, 3 F°, 1 D°, 3 F° excited 
states of carbon, carbon dimer, and naphthalene demon- 
strate the utility of our contribution. 

First, we reoptimized the primitive basis function ex- 
ponents for each system because the exponents of stan- 
dard QC and QMC basis sets, such as the Burkatzki, 
Filippi and Dolg (BFD) basis, represent a compromise. 



These standard exponents are designed to yield good en- 
ergies for some range of chemical environments and ex- 
citation levels, but they cannot be optimal for all sys- 
tems. We have shown that reoptimizing primitive basis 
function exponents for each system yields significant im- 
provements in the energy and fluctuations of the local 
energy, a. The most pronounced benefits were observed 
in higher-lying excited state calculations. In the most ex- 
treme case of 3 F° at the variational Monte Carlo (VMC) 
level, the 2z mixed basis was 30 mH lower in energy and 
110 mH lower in a than the 5z numerical basis. Al- 
though not discussed in this paper, we have found that 
reoptimization of standard Slater basis exponents used 
in all-electron calculations also provides considerable im- 
provements in energy and a. 

Second, we introduced Gauss-Slater (GS) basis func- 
tions for non-divergent pseudopotential calculations. GS 
functions behave like Gaussians at short distances and 
Slaters at long distances. In all systems considered, re- 
sults obtained using a mixed basis comprised of con- 
tracted and primitive basis functions improved when op- 
timized Gaussian primitives were replaced by optimized 
GS primitives. Importantly, for carbon dimer at the 
DMC level the 2z GS total energies are nearly converged 
with respect to basis size. 

A 3z mixed basis with optimized GSs for carbon atom 
or carbon dimer produces results comparable to the 5z 
BFD basis. Since the number of orbital coefficients to 
be optimized scales quadratically with basis size, the use 
of a more compact basis allows larger problems to be 
attacked in QMC. 
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Appendix A: Scaling Relations 

To derive the scaling relation between and , con- 
sider 

1 = (iV<) 2 / dr r 2 r 2 ^e- 2 ^y (Al) 



where 



c 



1 T (AK) 2 / duu 2 ^ 2 ^ (A2) 



2n+ 



1 fNS. 



(A3) 



C 2n+1 \K. 

Hence, the scaling relation for the normalization factor 
is 

N C = C+ l/2 N l_ (A4) 

The values of are given in Table IVII 

TABLE VI: Normalization factors for Gauss-Slater basis func- 
tions with unit exponent and principal quantum number n. 



n 


A 1 


1 


1.126467421 


2 


0.576609950 


3 


0.196581141 


4 


0.050275655 


5 


0.010280772 



To derive the scaling relations for the parameters a\ 
and cf in the Gaussian expansion of the Gauss-Slater 
functions, suppose the best-fit expansion for £ = 1 is 



Using Eqn. (|A4j) and performing the substitution r — > Qr 
results in 



ArC e i +<r 



~ C V V r(n + |) 



r(n + |) 



'2(2a£)"+i _c r = 



T(n+i) 



(A8) 



2„ 1 



Appendix B: Spatial Derivatives 



(A9) 
(A10) 



A general unnormalized radial basis function has the 
form 

R< n (r) =r n - 1 e g<( - r \ (Bl) 
where g^(r) is an arbitrary function. The gradient is 



dr 



where 



dr 

The Laplacian is 



(n-l) | dg<(r) 
r dr 



v 2 R c (r) d 2 Rj{r) 2dRj(r) 
n[ ' dr 2 + r dr ' 



where 



<9r 2 



+ 



aV(r) (n-i) 

dr 2 r 2 



1 /0i$(r) 



i£(r) V dr 
For Gauss-Slater functions, 

o , (CO 2 

BgC(r) r( 2 (2 + (r) 



dr 
d 2 g ( (r) 
dr 2 

For Gaussian functions, 



(1 + Cr) 5 
2C 2 

"(i + Cr) 3 ' 



/(r) = -Cr 2 
<9#C(r) 



dr 
d 2 g<(r) 
dr 2 



2(r 
2C- 



(A6) For Slater functions, 



dg^r) 



dr 

d 2 9 Hr) 
dr 2 



= 0. 



(B2) 



(B3) 



(B4) 



(B5) 

(B6) 
(B7) 
(B8) 

(B9) 
(B10) 

(BH) 

(B12) 
(B13) 

(B14) 
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Appendix C: Parameter Derivatives 



Appendix D: Exponents 



Wavcfunction optimization via the linear method re- 
quires both the derivatives of the wavefunction with re- 
spect to the exponent parameters £, and the Hamiltonian 
acting on those derivatives. From Eqn. (|A4I) . the deriva- 
tive of the normalization with respect to the exponent 
is 

«s . (j^±m Ni (C1) 

Now consider a general unnormalized radial basis func- 
tion of the form 



(C2) 



where g^(r) is an arbitrary function. The derivative of 
the radial part of the wavefunction with respect to the 
exponent is 

dRi(r) 



where 



The gradient is 
-dRi(ry 



d( 

fHr) 



fHrM(r), 



(C3) 



(C4) 



d( 

The Laplacian is 



df<(r) 
dr 



Ri(r)f + fHr)[VRi(r)] . (C5) 



2df<(r) 



r dr 
dr 2 

For Gauss-Slater functions 



dRjjr 
dr 



Ri(r) + f<(r)V 2 Ri(r). (C6) 



f C (r) 



(r 2 (2 + (r) 



(1 + CO 2 
d/ c (r) _ (r[A + Cr(3 + (r)} 

dr (1 + Cr) 3 

d 2 fi(r) _ 2C((r-2) 
dr 2 (1 + Cr) 4 ' 

For Gaussian functions, 

f C (r) = -r 2 
df<(r 



For Slater functions, 



dr 
d 2 fHr) 
dr 2 

f(r) 
dftjr) 
dr 
d 2 Hr) 
dr 2 



-2r 
-2. 



-1 
0. 



(C7) 
(C8) 
(C9) 

(CIO) 
(Cll) 

(C12) 

(C13) 
(C14) 

(C15) 



To promote use of this basis, exponents for each sys- 
tem studied are given in this appendix. Only exponents 
of the primitives are given, as the contracted functions 
are presented elsewhere [11]. Exponents for the ground 
state of carbon using CAS(4,13) CSFs, and, the ground 
state of carbon dimer using a single CSF are shown in Ta- 
ble IVIII Exponents for the lowest lying excited states of 



TABLE VII: Basis exponents for CAS(4,13) ground state of 
carbon and ground state of carbon dimer using G and GS 
basis functions. For each n, the nz basis includes n — 1 S 
primitives, n-lP primitives, n — 1 D primitives, and n — 2 
F primitives (where applicable). 



Type 


Size 


L 


C Exp. 


C 2 Exp. 


G 


2z 


S 


0.087 


0.145 






P 


0.129 


0.196 






D 


0.470 


0.679 




3z 


S 


0.102 


0.127 






s 


0.676 


0.998 






p 


0.104 


0.121 






p 


0.270 


0.423 






D 


0.314 


0.386 






D 


0.982 


1.099 






F 


N/A 


0.783 


GS 


2z 


S 


0.586 


0.853 






P 


0.984 


1.162 






D 


1.810 


3.774 




3z 


S 


1.000 


1.127 






s 


1.258 


1.570 






p 


1.059 


0.703 






p 


1.750 


1.416 






D 


1.132 


2.225 






D 


1.981 


3.228 






F 


N/A 


2.588 



carbon with 5 S°, 3 P°, 1 D° , 3 F° symmetries are shown in 
Table I Villi Exponents for the ground state of naphtha- 
lene are shown in Table ILXI In naphthalene, since atoms 
of the same atomic species are located at inequivalent 
geometrical locations, one could independently optimize 
the exponents for each inequivalent atom, but we have 
not done so because we expect the resulting gain to be 
small. 

We have found that the carbon S and P exponents 
change relatively little from one molecule to another 
(though they do differ more for the atom) while there 
is considerable leeway in the D exponents (they change 
considerably even from one optimization to another for 
a given molecule). This is because the energy and a are 
not as sensitive to the D basis functions as they do not 
appear in the ground-state determinant of the carbon 
atom. Hence it is possible to find an approximately op- 
timal set of exponents for the atoms in a large molecule, 
by optimizing them for a small molecule with the same 
atoms. 
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TABLE VIII: Basis exponents for the lowest lying excited 
states of carbon with 5 S°, 3 P°, 1 D°, 3 F° symmetries. For 
each n, the nz basis includes n — 1 S primitives, n — 1 P 
primitives, and n — 1 D primitives (where applicable). 



Type 


Size 


L 


'°S° Exp. 


S P° Exp. 


1 D° Exp. 


A F° Exp. 


G 


2 2 


S 


0.110 


0.006 


0.087 


0.112 






P 


0.150 


0.157 


0.109 


0.766 






D 


N/A 


N/A 


0.006 


0.006 




3z 


S 


0.356 


0.010 


0.094 


0.096 






s 


2.145 


0.284 


0.496 


0.689 






p 


0.809 


4.868 


0.079 


0.117 






p 


2.262 


6.707 


0.854 


0.524 






D 


N/A 


N/A 


0.007 


0.007 






D 


AT / A 

IN /A 


JN /A 


0.407 


0.822 


GS 


2z 


S 


1.960 


0.059 


0.756 


0.570 






P 


1.860 


3.002 


0.763 


0.588 






D 


N/A 


N/A 


0.351 


0.293 




3z 


S 


1.148 


0.371 


0.160 


0.894 






s 


1.451 


0.376 


0.931 


2.045 






p 


1.137 


0.691 


0.043 


1.326 






p 


1.325 


2.087 


0.444 


3.016 






D 


N/A 


N/A 


0.238 


0.295 






D 


N/A 


N/A 


0.953 


1.545 



TABLE IX: Basis exponents for ground state of naphthalene, 
CioHg, using G and GS basis functions. For carbon, the 2z 
basis includes 1 S primitive, 1 P primitive, 1 D primitive. 
For hydrogen, the 2z basis includes 1 S primitive. 



Type 


Size 


L 


CioH 8 C Exp. 


CioH 8 H Exp. 


G 


2z 


S 


0.139 


0.099 






P 


0.191 


N/A 






D 


0.754 


N/A 


GS 


2z 


S 


0.875 


0.798 






P 


1.118 


N/A 






D 


2.109 


N/A 



